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Abstract 

Transition from chaotic to quasi-periodic phase in modified Lorenz model is analyzed 
by performing the contact transformation such that the trajectory in i? 3 is projected 
on R 2 . The relative torsion number and the characteristics of the template are 
measured using the eigenvector of the Jacobian instead of vectors on moving frame 
along the closed trajectory. 

Application to the circulation of a fluid in a convection loop and oscillation of 
the electric field in single-mode laser system are performed. The time series of the 
eigenvalues of the Jacobian and the scatter plot of the trajectory in the transformed 
coordinate plane X — Z in the former and \X\ — \Z\ in the latter, allow to visualize 
characteristic pattern change at the transition from quasi-periodic to chaotic. In the 
case of single mode laser, we observe the correlation between the critical movement 
of the eigenvalues of the Jacobian in the complex plane and intermittency. 
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1 Introduction 



In dissipative dynamical systems exhibiting transition from periodic or quasi- 
periodic behaviour to chaotic one, topological method is a useful tool. Gilmorefl] 
has recently reviewed how the integer invariants can be obtained from chaotic 
time series. In the time series, sudden bursting is called intermittency. Pomeau 
and Manneville[2,3] proposed three types of routes to chaos which accompany 
intermittency. The type I is due to disappearance of a stable periodic orbit 
via pitchfork bifurcation, the type II is via period doubling bifurcation and 
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the type III is via Poincare-Andronov-Hopf bifurcation. The definitions of the 
type of bifurcation can be seen in e.g. [11,19]. 

In the Lorenz model, which simulates the flow of fluids between hot lower 
plate and cold upper plate, chaotic behaviour arises at a certain parameter 
region through the pitchfork bifurcation. Birman and Williams [4] showed that 
the strange attractor of the Lorenz model in R 3 can be constructed on a 2- 
dimensional bifurcating manifold which is called knot-holder or a template, 
and that the topological analysis of the orbits can be reduced to that of braid 
groups. 

By applying external temperature oscillation p\ cos out on the upper plate, 
one can stabilize the chaotic behaviour in the Lorenz model. In this case the 
dynamics is extended from R 3 to R 3 x S* 1 where S* 1 is the 1-dimensional 
torus, and one can guide the trajectory of the orbit from chaotic to periodic 
or quasi-periodic one. 

Holmes and Williams [5] extended the Lorenz mapping to mappings that con- 
tain Smale's horseshoe mapping which is characterized by stretching and fold- 
ing, and expressed dynamics on the horseshoe template as a sequence of sym- 
bols L and R corresponding to the relative position of the two orbits on the 
template. They showed that the period-doubling bifurcations can be speci- 
fied by the number of half twist £ and the index of folding p and called the 
associated 2-dimensional manifold as (£,p)-template. 

Aizawa and Uezu[6-8] performed the topological analysis of the perturbed 
Lorenz model and specified the bifurcating orbits by the linking number of 
two orbits L(C±, C2), torsion number rii and relative torsion number r\ which 
is equivalent to p/2 of the Holmes- Williams, and derived a recursion formula 
for those of the 2 fc th bifurcation orbit rf^ and n\ k K 

Arimitsu and Motoike[9-ll] extended the analysis of Holmes and Williams 
and showed that from the power spectrum of the bifurcating orbits, topolog- 
ical properties of the orbit can be extracted. They defined the local crossing 
number C/ n \, global crossing number ci n \ and the linking number l/ n ,n-i) of the 
2™th and the 2 n_1 the bifurcation orbit and obtained the recursion formula. 

The chaos control of the perturbed Lorenz system was studied from the tech- 
nological point of view[12,13]. By a simple picosoot type parametric driving 
of a convection loop model, it was shown that the chaotic system after driven 
to periodicity can be led back to chaotic behaviour through intermittency. 
General methods of chaos control are proposed by several authors[14,15]. 

There are analogies between the Lorenz model and the single mode laser sys- 
tem[16,17]. The dynamics of the electric field, polarization and the inversion 
density can be assigned as the coordinate of the three dimensional manifold, 
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but the electric field and polarization are accompanied by the modulation, 
and they are expressed as complex numbers. The transition to chaos by ap- 
plying external constant electric field to the single mode laser system was 
simulated [18]. 

In this paper, we analyze the convection loop model and the single mode 
laser system from a slightly different point of view. We adopt the framework 
of differential geometry and adopt the contact transformation to simplify the 
analysis. We consider the topological properties of these models following [7, 10] 
and consider an extension of our method to other systems. 



2 Contact transformation 



When an orbit in a d-dimensional dynamical system is specified byt,xi,X2,---, Xd, 
Pi,P2, • • • iPd-, and if the coordinate transformation 

T = T(t,x u x 2 , - ■ ■ ,x d ,p 1: p 2 , - • • ,p d ) (1) 

Xj = X j (t,x 1 ,x 2 , - ■ • ,Xd,Pi,P2, - • ■ ,Pd) (2) 

Pj = P j (t,x 1 ,x 2 ,-'-,Xd,Pi,P2,---,Pd) (3) 

does not change the total differential equation 

dt — p\dx\ — p 2 dx 2 — • • • — p d dxd = (4) 

i.e. when the equation 

dT - P 1 dX 1 - P 2 dX 2 P d dX d 

= p(t, x u x 2 , ■ ■ ■ , x d , pi, p 2 , • • ■ , p d ) {dt - pidxi - p 2 dx 2 Pddx d ) 

= (5) 

is identically satisfied, the transformation is called contact transformation. 

When the dynamical equation is given as 

^ = h 1 (x,y) 

— = h 2 {x,y,z) 
dz 

— =h 3 (x,y,z) (6) 
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dx 

we perform the transformation X = x, Y = hi(x,y). By using — = hi(x,y) 

dy dY 
and — = h 2 (x,y,z), the derivative can be expressed as g 2 (X,Y, z), and 

we define it as Z. Further using — = h 3 (x,y,z), the derivative — — can be 

dt dt 
expressed as f(X, Y, Z). Thus the system can be expressed as 



= Y 



dX 
~dt 

dt 



(?) 



To show that it is a particular case of contact transformation, it is sufficient 
to consider the identity 

Zf(X, Y, Z)dX + Yf(X, Y, Z)dY + YZdZ 

YZf(X,Y,Z) 
h 1 (x,y)h 2 (x,y,z)h 3 (x,y,z) 

x (h 2 (x, y, z)h 3 (x, y, z)dx + h^x, y)h 3 (x, y, z)dy + h ± (x, y)h 2 (x, y, z)dz) 
= p(x, y, z)(h 2 (x, y, z)h 3 (x, y, z)dx + hi(x, y)h 3 (x, y, z)dy 
+h 1 (x,y)h 2 (x,y,z)dz) (8) 

When the system is discretized we obtain 



X n+ i = X n + Y n dt 
Y n+ \ = Y n + Z n dt 
Z n +i — Z n + f(X n , Y n , Z n )dt 



(9) 



The Jacobian of this mapping is 



/ d(X n+1 -X n ) d(X n+1 -X n ) d(X n+1 -X n ) \ 

dX n dY„ dZ n 

d(Y n+1 -Y n ) d(Y n+1 -Y n ) d(Y n+1 -Y n ) 

8X n dY n dZ n 

d(Z n+ i — Z n ) d(Z n+ i — Z n ) d(Z n+ i—Z n ) 

\ dX n dY n dZ n 

( 



/dt 








1 







1 



df(X n+ l,Yn + l,Z„+l) df(X n + l,Yn+l,Z n + l) 8f(X n+ l ,Yn + l ,Z„ + i) 

8X n 8Y n dZ n 



(10) 
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and its eigenvalues are given by the solution of 



\3 9f(X n+ i, Y n+ i, Z n+ i) 2 df(X n+1 , Y n+ i, Z n+ i) 



df(X n+ i, Y n+1 , Z n+ i) _ q 



dZ n 

which means that there is a zero eigenvalue if — — — — VdlL = o is 

oZ n 

satisfied. 



3 Lorenz equation 

In the case of Lorenz equation 



dx 
dy 

— = —y + rx — xz 
dt 

dz 

— = -bz + xy (12) 



dZ 

The derivative — = f(X, Y, Z) becomes 

(Jjv 



f(X, Y, Z) = (YZ + aY 2 + Y 2 - aXZ -XZ - X 3 Y - aX 4 

-bXZ - abXY + abrX 2 - bXY - abX 2 )/X, (13) 



3.1 Scatter plot in the X — Z plane 



We consider a convection loop of water which has the Plandtl number a = 
10 and the geometrical factor b = 8/3. The threshold r t for these values is 
24.7368[19]. The time series of x for r t = 24 is shown in Figure 1 and that for 
r t = 25 is shown in Figure 2. We find clearly that the system converges to a 
fixed point when r < r t but intermittency occurs when r > r t . 

After the contact transformation the scatter plot in the X — Z plane is shown 
in Figures 3 and 4. 
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We remark that in the X-Z plane, the scatter plot X(Z) changes from the 
monotonic function in the case of r < r t to 3-valued function in the case of 
r > r t . 



4 Convection loop model 



We consider the stability of thermosyphon or circulation of water in convection 
loop with heat source at the bottom and cooler at the top[13]. At the cooler, 
modulation of temperature in the form p + Pi cos wt is assumed. The modified 
Lorenz equation is 



^ = a(y-x) (14) 
dy 

— = —y + (po + Pi cosutjx — xz (15) 
at 

dz 

— = —bz + xy (16) 
at 

Due to pi term, quasi-stable orbit can appear even when p > r t . 
dZ 

The derivative — = f(X, Y, Z) after the contact transformation is 



/(X, Y, Z, t) = (YZ + aY 2 + Y 2 - aXZ - XZ - X 3 Y - aX 4 

-bXZ - abXY + ab(p + pi coscot)X 2 - bXY - abX 2 )/X 
— cr(piu) sin u>t)X, (17) 

4-1 Scatter plot in the R 3 and in the X — Z plane 



In Figures 5 and 6, we show quasi-periodic flow in R 3 for p = 26, p\ = 
2.5, uj = 9 and for p = 28, pi = 11, cu — 8.5, respectively. After the contact 
transformation, the corresponding quasi-periodic flows in the X-Z plane are 
shown in Figures 7 and 8, respectively. 

We assign eigenvalues of the Jacobian (10) as Aj, i=l,2,3 such that ReAi < 
ReA 2 < ReA 3 . The branching of the eigenvalues in the quasi-periodic region are 
shown in Figure 9(a), 10(a), for the small pi = 2.5 system and large p 1 = 11 
system, respectively. We observe that, when the steps are started from near 
the origin, three real roots appear in the beginning but subsequently complex 
conjugate pairs corresponding to rotation around two fixed points dominate. 
In the plot of time-series of eigenvalues, there appears a discontinuity when 
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X = 0. They are absent when contact transformation is not performed and 
do not affect the measurement of relative torsion number. 



4-2 Topological analysis 

For a dynamical system defined by 

^ = F{x,t,n) (18) 

where \x is the bifurcation parameter, one defines eigenvalues of 
OF 





r OF 

S(t,0) = exp[J — (x (t'),n)dt'} (19) 



as Xi, (i = 1, 3), such that@|Ai] < | A2 j < | A3 1 , and eigenvectors associated with 
them as (i — 1, 3). Around the periodic orbit 

C = {x (t)\0<t<T}, (20) 



one defines two closed orbits [7] 

C K = {x (t) + Cwi(t)\0 < t < mT} (21) 

where ( is a sufficiently small number and Wi is the normalized eigenvector 
w l (t) = e, l (t)/\\e l (t)\\. (22) 

The parameter m is 1 when A« > and 2 when \ < 0. 

In the system after the contact transformation, we define the normalized eigen- 
vectors of Jacobian where X = (X, Y, Z), but not that of S(t, 0) as Wi. 

The local torsion number r(w), which provides information on the rotation 

number of Ci,(=i around Co is defined by using f 1 = — : — , f 2 = — = — , 

\\X || || X x X || 

/ 3 = fi x /2 and a j(t) =w(t)- fj as 

r(w) = J- / a2it) ' - ' a3(t) (23) 



2nJ a 2 (t) 2 + a 3 (t) 
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In Figure 9(b,c), and 10(b,c) we show the third component of the two nor- 
malized eigenvectors Wi and w 2 for the small pi system and large pi system, 
respectively. The horizontal lines of the third component of the eigenvectors 
Wi and w 2 equal 1 indicate stability of the orbit C . In the case of small pi 
(Figure 9(b)), we find two stable windows, which indicates LR type period 2 
structure. While in the case of large pi (Figure 10(b)), there is an additional 
breaking in the eigenvector W\ which corresponds to a compensating twist on 
the LR type period 2 structure. 

The eigenvector w 3 has also small first and second components, and its third 
component oscillates between -1 and +1. We measure the relative torsion 
number of the orbit C^i, since Wi and w 2 specify the trajectory of Co- We 
observed that the local torsion number r(w 3 ) in the case of small pi is -1/2 
and that of large pi is 1/2. 

In order to get local crossing number Cj of the orbit [10], we measure the 
position of the highest peak in the power spectrum of the orbit of Figure7 and 
Figure8. The results of small pi and large p\ together with the pure chaotic 
case without perturbation (pi = 0) are shown in Figures 11,12 and 13,14 
respectively. Templates are characterized in terms of number of half twists 
£ and the direction of folding p e {—1, 1}, as (£,p), and the local crossing 
number C( 2 ) = 2£ + p in the case of period doubling bifurcation [10]. 

In the power spectrum of po = 26 and pi = 2.5, ui = 9, we observe the peak 
position is at 00 = 4 in arbitrary unit. When the perturbation is set to pi = 0, 
the position moves to u — 3. Hence, using the information r(w 3 ) = —1/2, 

C 2 3 

i.e. p = —1, we assign — = -. Since C 2 = 3, the template of small p 1 is 
&P) = (2,-1). 2 

In the power spectrum of p = 28 and pi = 11, uj 2 — 8.5, the peak position 
at uj = 4 and it does not change after setting p 1 = 0. Using the information 

r(w 3 ) = 1/2, i.e. p = 1, we assign C\ = -. The template of large pi is 

&P) = (0,1). 



5 Single mode laser 

In 1975, Haken showed that under certain conditions a coherently pumped ho- 
mogeneously broadened ring laser obeys the Lorenz equation[16,17]. Whether 
laser actually behaves in this way was studied, and data of a kind of infrared 
single-mode gas laser was simulated [18]. 

In this model, the inversion density is defined as D, the electric field and 
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polarization are specified as x = xe ldKt and p = pe tAt , respectively. In the 
single mode laser, A = 9k, and the coupled equation of the system with 
external static electric field y becomes 

dx 

— = -kiiOx + x-y + 2Cp) (24) 
dt 

^ = xD-(l + iA)p (25) 
^ = -l\\{xp* +x*p) + D + l] (26) 



We consider the modulation g = e lAi as a structure group U(l) defined on the 
base space. The complex space specified by x,p and D is projected onto -R 3 
specified by \x\, \p\ and D by identifying all the fields x and p that are related 
as g~ x xg = x'. 

The stable solution that is obtained by taking l.h.s. of the equations to be zero 
gives the relation[18] 

^i[(i- TT #^) 2 +c+ TT f^)r 2 (27) 



The function \x\(y) shows a bending behaviour and in the region 300 < \x\ < 
500 it is 3- valued function(Fig.l5). 

On the projected space we can perform the contact transformation as before 
and we obtain 



dZ 

—— = —inOZ — kZ 
dt 

-2kC[YF l + X(-7(^(XF 2 * + X*F 2 ) + F 1 + 1))] 

-(Z + ik6Y + kY)(l+iA) (30) 

where 

-Z-ik6Y-kY-(Y + ikex + kX-ky)(l + iA) , 
Fl = 2kCX (31) 
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5.1 Scatter plot in the \X\ — \Z\ plane 

Due to the positive real y, the system is unstabilized and shows a strong 
sensitivity to the initial condition. In this paper, we present the result of 
k = \,C = 500, A = 9 = 5, and the initial condition X = 350, Y = 0, Z = 1. 
The scatter plot of \Z\ in the \X\ — \Z\ plane for each y from 230 to 340 are 
plotted in Figure 17 for y > 270 and Figure 16 for y < 260. Convergence to 
a stable branch can be observed for y > 300, but at y — 270 intermittency 
occurs and the situation becomes chaotic for y < 260. 

Although the chaotic behaviour in the time series of x can be observed by the 
bursting and spiking phenomena [18], the scatter plot in the \X\ — \Z\ plane 
shows clearer transition between chaos and quasi-stable mode. We observe 
intermittency begins to occur at around y = 270 and the system become 
chaotic for smaller y values. 

• Fig. 16 Scatter plot of (|-X"|, \Z\) for y < 260. The points correspond to 
y=230, 240, 250, 260, from left to right. See the extra file: laserl.gif 

• Fig. 17 Scatter plot of {\X\, \Z\) for y > 270. The points correspond to 
y=270, 280, • • •, 340, from left to right. See the extra file: laser.gif 

5.2 Intermittency 

The scatter plot of (|-X"|, \Z\) at y — 270 is presented in Figure 18. The inter- 
mittency near \X\ = 53 is shown in Figure 19. 

In order to check the type of intermittency, we measure the time series of 
the eigenvalues of the Jacobian. Figures 20, 21 and 22 are the scatter plot 
of the time series of the three complex eigenvalues. There appears a pair of 
boomerang type and a ring type. 

As the point (|-X"|, \Z\) moves from right, makes a loop, and then move to left 
as shown in Figure 19, the eigenvalues of ring type move from around (0, —3i) 
counter clock wise, move to around (—0.5, — 6.5i), make a cusp, and then move 
clock wise as shown in Figure 22. The movements of eigenvalues of boomerang 
type are counter clockwise. 

As the eigenvalue of boomerang type approach the real axis, \Z\ takes a local 
maxima, and as the eigenvalue of ring type makes a cusp, \Z\ takes a local mini- 
mum. These properties can be qualitatively understood as the movement along 
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the trajectory of the U(l) group makes a jump of the flow in (\X\, \Y\, \Z\) base 
space, which has the multivalued structure as in the original Lorenz model. 



6 Discussion and outlook 

We applied the contact transformation to the Lorenz equation and its exten- 
sion, and analyzed transition from chaos to quasi-periodic phase in convection 
loop and in the single mode laser. The branching of eigenvalues of the Jaco- 
bian after contact transformation allows visualizing transition from orbit near 
one fixed point to the other one. 

In the case of convection loop we considered perturbation by p\ cos out term 
with large p\ and small p\. From the behaviour of eigenvectors of the Jacobian, 
we observe transition from quasi-stable to unstable orbits. The relative torsion 
number can be measured by using the eigenvectors and the local crossing num- 
ber can be measured from the power spectra of Z(t) in periodic orbit with 
pi cos cut type perturbation and without perturbation. The topological data 
after contact transformation are simpler than those before contact transfor- 
mation. Without contact transformation, the measurement of relative torsion 
number is relatively ambiguous since for x, y and z directions are equiv- 
alent and there is no guideline for selecting one like Z. 

The X — Z plot in the convection loop and \X\ — \Z\ plot in the single mode 
laser allows to visualize characteristic pattern change when the system changes 
from chaotic to quasi-periodic phase. 

We observed intermittency in the single mode laser system at around y = 270 
and observed correlation between the critical movement of complex eigenvalues 
of the Jacobian and the sudden large \Z\ oscillation. The orbit in the fiber 
space specified by the U{1) group causes a jump in the trajectory in the base 
space specified by (|X|, \ Y\, \Z\). It is a new type of intermittency. 

The topological characteristics of the orbit can be attributed to the jump in 
the relative linking number of the orbit Xq+w 3 around the quasi-periodic orbit 
Xq, where w 3 is the normalized eigenvector of the Jacobian, associated with 
the eigenvalue possessing the largest real part. We observe that the behaviour 
of the eigenvectors corresponding to eigenvalues with the smallest and the 
second smallest real part give information of the template and the orbit Co- 

When the Lyapnov exponents satisfy Ai > A 2 = > A 3 > • • • > A^ and there 
is a strong local attraction, the Lyapnov dimension becomes 

d L = 2 + A- < 3 (33) 
I A3 1 
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which means that the system can be mapped on 2-dimensional template [20, 21]. 
In the decomposition of the three dimensional space to two dimension + one 
dimension space, the 2-dimensional template could always contain stable man- 
ifold. 

In higher dimensional space, 4-dimension as an example, the chosen 2-dimensional 
template could contain only the unstable manifold. To study it, we consider 
the extended Lorenz equation[22] in R 4 . 



dx 

^ =*(*-*) 

dy , \ 

— = — y + rx — x{z + w) 
dt 

dz 

Tt =-bz + xy 

d ^ = kx 2 -cw (34) 
dt 

where parameters are a — 16, b — 4, c = 1, k — 0.02. 



r - 1 

In this model, fixed points exist at x — y — ±w — and a pitchfork bifur- 

V k + l/o 

cation takes place for the parameter r ~ 530 [22]. Around this point kx 2 — cw 
becomes positive and w increases monotonically but in the projected (x, y, z) 
space we can observe the attractor. 

When we perform the contact transformation 

« y, £ =z, f = W (35) 
dt dt dt v 1 



— = (Z(Y + Zf - a(WX 2 ((l + b + c)X - 2Y) + cX 5 Y - cX 2 Y 2 + X A Y 2 
dt 

+ 2XY 3 + cX 3 Z + X 5 Z - 2X 2 YZ - cX 2 YZ - 2Y 2 Z 

+ 2XY 2 Z - 2YZ 2 + bX 2 (-Y 2 + XZ-YZ + cX(Y + Z)) 

- a 2 (WX 3 + bkX 6 + 2X 5 Y + 2kX 5 Y - bX 2 Y 2 + 2XY 3 + bX 3 Z - 2X 2 YZ 

- Y 2 Z + cX 2 ((-l + r)X 2 + X 4 - Y 2 + X(bY + Z)))/(aX 3 ) (36) 

we find spiral divergence in the Z — W plane. 

The dynamical flows in these frames are trapped in their corresponding un- 
stable manifolds. To remedy it, one can choose w as X and linearize kx 2 — cw 
as k(x^ + Xqx) — cw where x is the coordinate of the fixed point and study 
dynamics on the template in the co-dimensional space of w, i.e. x, y, z space. 
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In general, when the dimension of the space is larger than three, it is neces- 
sary to choose coordinate in the contact transformation so that the template 
contains the 2-dimensional stable manifold. 

I acknowledge with thanks that graduate students M. Horinouchi and N. 
Kamimura obtained some simulation data in this paper. Thanks are also due 
to the referee for bringing references [6-8, 22] to my attention and Dr. T. Arim- 
itsu for kindly sending me his references and helpful discussion. 
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Fig. 1. Time series of x. r = 24, a = 10, 
6 = 8/3. 
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Fig. 3. Scatter plot in X-Z plane, 
r = 24, a = 10, 6 = 8/3. 




Fig. 5. Perturbed Lorenz model. 
po = 26, pi = 2.5, oj = 9. 
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Fig. 7. Scatter plot in X — Z plane. 
po = 26, pi = 2.5, uj = 9, 
4750At < t < 5850At, At = 0.0015 
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Fig. 2. Time series of x. r = 25, a = 10, 
6 = 8/3. 



« 




• 





Fig. 4. Scatter plot in X-Z plane, 
r = 25, a = 10, 6 = 8/3. 




Fig. 6. Perturbed Lorenz model. 
p = 28, pi = 11, lo = 8.5. 
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Fig. 8. Scatter plot in X — Z plane. 
po = 28, pi = 11, u = 8.5, 
4000At < t < 5600At, At = 0.001 
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Fig. 9. Time series of the real part of the eigenvalues of Jacobian(a), the third 
component of the eigenvector corresponding to the lowest eigenvalue of Jacobian(b) 
and that corresponding to the second lowest eigenvalue, po = 26, p\ = 2.5, b = 8/3, 
u = 9. 
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Fig. 10. Time series of the real part of the eigenvalues of Jacobian(a), third com- 
ponent of the eigenvector corresponding to the lowest eigenvalue of Jacobian(b), 
and that corresponding to the second lowest eigenvalue, r = 28, p\ = 11, b = 8/3, 
u = 8.5. 
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Fig. 11. The power spectrum of Fig. 12. The power spectrum of 
po = 26, pi = 2.5, lo = 9. p = 26, pi = 0. 
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Fig. 13. The power spectrum of Fig. 14. The power spectrum of 
po = 28, pi = 11, u = 8.5. po = 28, pi = 0. 
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Fig. 15. The electric field \x\ of the sta- 
ble or quasi stable orbit function 
of the external static field y. 
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Fig. 19. Scatter plot of m 
so ioo 150 200 250 300 350 the case of y = 270 near \X\ = 53. 

The right local maximum, local mini- 
mum and the left local maximum are 
assigned by circles. 



Fig. 18. Scatter plot of (\X\, \Z\) in the 
case of y = 270. 
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Fig. 20. Scatter plot of complex eigen- 
values of the Jacobian. Eigenvalues of 
boomerang type which have negative 
real part. Ordinate is the real part 
and abscissa is the imaginary part. The 
three circles from right to left of Figure 
19 correspond to the circle in the right, 
in the top and in the left, respectively. 
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Fig. 21. Scatter plot of complex eigen- 
values of the Jacobian. Eigenvalues of 
boomerang type which have positive 
real part. The three circles from right 
to left of Figure 19 correspond to the 
circle in the left, in the bottom and in 
the right, respectively. 
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Fig. 22. Scatter plot of ring type com- 
plex eigenvalues of the Jacobian. The 
three circles from right to left of Figure 
19 correspond to the circle in the 2nd 
quadrant, at the cusp and in the 3rd 
quadrant, respectively. 
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